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Abstract. We describe the anisotropy of dark matter clump distribution caused by 
tidal destruction of clumps in the Galactic disk. A tidal destruction of clumps with 
orbit planes near the disk plane occurs more efficiently as compared with destruction 
of clumps at near-polar orbits. A corresponding annihilation of dark matter particles 
in small-scale clumps produces the anisotropic gamma-ray signal with respect to the 
Galactic disk. This anisotropy is rather small, ~ 9%, and superimposed on that due to 
off-centering position of the Sun in the Galaxy. The anisotropy of annihilation signal 
with respect to the Galactic disk provides a possibility to discriminate dark matter 
annihilation from the diffuse gamma-ray backgrounds of other origin. 
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1. Introduction 

A primordial power-law spectrum of density fluctuations in the Dark Matter (DM) 
ranges from the largest scales above the scales of superclusters of galaxies to the smallest 
sub-stellar scales according to prediction of inflation models. This permits to predict 
the properties of smallest DM structures from the known CMB fluctuations at large 
scales. Substructures of DM in the galactic haloes with a rather large mass, > 1O 7 M , 
were extensively discussed in early works, see for example pQ. The nonlinear dynamics 
and mechanism of hierarchical clustering of these large DM clumps were analyzed in 
both analytical calculations [21 [3, H] and numerical simulations [5j El [7] . At sub-stellar 
mass-scales of DM fluctuations, a principal new phenomenon arise — the cutoff of 
mass spectrum due to collisional and collisionless (free streaming) damping processes 
of DM particles in the forming clumps. The resulting smallest mass of DM clumps is 
determined by the properties of DM particles, in particular, by their elastic scattering. 
For detailed calculations of this cutoff see e. g. [8] and references therein. Additionally 
the cutoff of mass spectrum is influenced by the acoustic absorption [9] at the time 
of kinetic decoupling of DM particles [TO] and also by the horizon-scale perturbation 
modes [11] . In [12] the kinetic equations for DM phase space density were solved in 
the case of perturbed cosmological background by taking into account the acoustic 
absorption, horizon-scale modes and gravitational perturbations. A corresponding value 
of the smallest clump mass for neutralino DM is of the order of the Moon or Earth mass. 

The formation of small-scale DM clumps with a mass larger than the Earth mass, 
M m i n ~ 10 _6 M Q , have been explored in numerical simulations [131 EE]- A resulting 
differential number density of small-scale clumps, n(M) dM oc dM/M 2 , turns out very 
close to that obtained in the numerical simulations of large-scale clumps with mass 
M > 1O 6 M . The other important result obtained in numerical simulations [13] is 
determination of the internal density profile in the isolated clump of minimal mass. The 
resulting density profile is approximately a power-law, oc r _/3 , with /3 = 1.5 — 2.0, which 
is in a good agreement with theoretically predicted value (3 = 1.7 — 1.8 according to [2]. 

The number density of small-scale DM clumps existing nowadays in the Galactic 
halo is determined by their tidal destruction during hierarchical structure formation 
[15] and also by tidal interactions with stars in the Galaxy [T3l [TBI [T71 [181 [T9l 1201 [21] . 
Annihilation of DM particles in small-scale clumps [HI [221 [231 (2H [251 ES [271 12HI [29] 
enhances the total DM annihilation signals in our Galaxy and thus boosts a chance for 
indirect detection of DM. 

The usual assumption in calculations of DM annihilation is a spherical symmetry 
of the Galactic halo. In this case an anisotropy of annihilation gamma-radiation is only 
due to off-center position of the Sun in the Galaxy. Nevertheless, a principal significance 
of the halo nonsphericity for the observed annihilation signal was demonstrated in [30J. 
According to observations, the axes of the Galactic halo ellipsoid differ most probably 
no more than 10 — 20%, but even a much more larger difference of axes, up to a factor 2, 
can not be excluded [311 [32]. This leads to more than an order of magnitude uncertainty 
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in the predicted annihilation flux from the Galactic anti-center direction [3D]. It must 
be noted also the "intrinsic" annihilation anisotropy caused by the small-scale DM 
clustering itself. A corresponding angular power spectrum of annihilation signal at 
small scales is connected with a power spectrum of DM clumping [33] . In principle, the 
DM clumps may be seen as point sources at the gamma-sky [30J. Another minor source 
of annihilation anisotropy is a dipole anisotropy due to proper motion of the Sun in the 
Galaxy [M]. 

In [IT] the anisotropy with respect to the Galactic disk was discussed basing on the 
numerical calculations of the destruction of DM clumps by stars in the disk and taking 
into account the influence of gravitational potential of the disk on the clump orbits. 
It was also shown [151 [19] that (i) small-scale DM clumps dominate in the generation 
of annihilation signal and (ii) the Galactic stellar disk provides the main contribution 
to the tidal destruction of clumps at r > 3 kpc, i. e. outside the central bulge region. 
A process of clump destruction in the halo is anisotropic in general (e. g. it depends 
on the inclination of clump orbit with respect to the disk plane). Respectively, the 
DM annihilation in the halo is also anisotropic. In this work we estimate the value of 
this anisotropy. It must be stressed that with a present state of art it is impossible to 
separate this source of anisotropy from that produced by the halo nonsphericity. More 
detailed investigation is required to constrain the shape of the halo and to search the 
distinctive features of annihilation anisotropy due to non-spherical halo clumpinesss. 
The detectors at the GLAST satellite will be sensitive to anisotropy up to 0.1% level 
[34] . This will provide a hope to discriminate the anisotropic DM annihilation signal 
from the diffuse gamma-ray backgrounds. 



2. Destruction of clumps by disk 

Crossing the Galactic disk, a DM clump can be tidally destructed by the collective 
gravitational field of stars in the disk. This phenomenon is similar to the destruction of 
globular clusters by the "tidal shocking" in the Galactic disk [35]. The corresponding 
energy gain per unit mass of a clump at one disk crossing [35] is 



AE = 2 ^\ (1) 

vi 



z,c 



where g m is the maximum gravitational acceleration of the clump moving through the 
disk, Az is a vertical (perpendicular to the disk plane) distance of a DM particle from the 
clump center, v ZjC is a vertical component of velocity at disk crossing. The dependence 
of v 2 jC on the inclination of orbit relative to the disk plane is the origin of the discussed 
anisotropy in the clump destruction, and, as a result, the origin of the anisotropy in 
annihilation signal. 

The surface mass of the Galactic disk [36] can be approximated as 
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with Md = 8 x 10 10 M Q and r = 4.5 kpc, and therefore 

g m (r) = 2nGa s (r). (3) 
We use the power-law parametrization [2j [3j IU [13] of the internal density of a clump 

Pint(r) = ^— p - , (4) 



3 r \R, 

where p and i? are the mean internal density and a radius of clump, respectively, /? = 1.8 
and /Oint(r) = at r > i?. The total (kinetic plus potential) internal energy of a clump 
for density profile (jlj) is given by 

iff, 3-/3 

~ 2(5 - 2/3) i2 ' (5) 
where M is the mass of the clump. Integrating (PQ) over a clump volume and using the 
density profile (T4|), one obtains an energy gain for the whole clump as 

AE (5 - 2/3) g 2 



> in 



\E\ tt(5- P)Gpvl c ' ^ 

We will use the following criterium for a tidal destruction of clump: a clump is destructed 
if a total energy gain AEi after several disk crossings exceeds the initial internal energy 
\E\ of a clump. 

Let us consider now some particular orbit of a clump in the halo with an 
"inclination" angle 7 between the normal vectors of the disk plane and orbit plane. 
The orbit angular velocity at a distance r from the Galactic center is dip/dt = J / (mr 2 ), 
where J is an orbital angular momentum of a clump. A vertical velocity of a clump 
crossing the disk is 

w 2 , c = ^-sin7, (7) 

mr c 

where r c is a distance of crossing point from the Galaxy center. There are two crossing 
points (with different values of r c ) during the one orbital period. The momentum 
approximation used here for calculations of the tidal heating is violated at small 
inclination angles, 7 -C 1. Nevertheless the resultant anisotropy is a cumulative quantity. 
It results from an integration over all clump orbits, and orbits with 7 1 provide only 
small input into the anisotropy value. 

We use here the Navarro- Frenk- White (NFW) density profile for the Galactic halo: 

PH(r) = (r/L)(i° + r/Lr (8) 
where po is normalized in such a way that at position of the Sun Ph(^o) = 0.3 GeV/cm 3 
and L = 45 kpc. 

Let us introduce the following set of dimensionless variables: 

L po Stt&PqL^M 2 

E orh /M - $0 $ - $0 

£ AttGpoL* ' P AnGpoL^ [ } 
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where $ = —^GpoL 2 , E orh is an total orbital energy (kinetic and potential) of a clump. 
In these dimensionless variables the density profile of the halo (jSJ) has a form 

Consider the Galactic halo model with an isotropic velocity distribution, which is 
appropriate for the halo formed by the hierarchial clustering of sub-haloes. In this 
model the distribution function depends only on energy. The gravitational potential 
ip(x), which corresponds to the density profile (iTTj) is given by 

1,(x) = 1 - lQg(1 + X) . (12) 

x 

The equation for the turning points of an orbit, r 2 = 0, for the potential ( 1T21) can be 
written as 

log(l + x) y . , 

1 - — + 4 = e - 13 

x x z 

From ffl3|) one can find numerically the minimum x min and maximum x max distance 
between a clump and the Galactic center as function of e and y. Now denoting p = cos 9, 
where 9 is an angle between the radius- vector r and the particle velocity v, we have 

y={l-p 2 )x 2 {e-^{x)). (14) 

According to our assumption, the vectors v/v are distributed isotropically at each point 
x, and therefore p has a uniform distribution in the interval [0, 1]. The relation between 
p(x) and the distribution function F(e) is given, according to [37], by 

i 

p(x)=4ttV2 J de[e - ip(x)} 1/2 F(e). (15) 
The function F(e) for a halo profile (iTTj) can be fitted [38] with a good accuracy as 



F(e) = F 1 (l - ef 2 e" 5 / 2 



ln(l-e) 



e 

1-2 



1 



(16) 



where Fi = 9.1968 x 1(T 2 , P = EPi(l - e)\ (q,Pi,P2,P3,Pi) = (-2.7419,0.3620, - 
0.5639,- 0.0859,- 0.4912). 

The time of motion from x m j n to x max and back is 

^ 3^ max , 

T c (g,e,p)= ; = / , 17 



where y depends on x, e and p according to (fl4"j) . The trajectory (the orbit) of a clump 
in the potential ffl2|) is not closed. A precession angle during the time T c /2 is 

* = yl/2 / 2/ 7H n _ (18) 

and < 0. Therefore the orbital period is greater then T c and equals to 

T t = T c (l+0/7r) _1 . (19) 



Anisotropy of DM annihilation 6 

During a life-time of the Galaxy, £q ~ 10 10 yrs, the relative total energy gain of a clump 
due to the tidal heating in multiple crossings of the disk is 

AF 1 N 

^ = ^E( A ^ + A ^), (20) 

where AE\ and AE % 2 are given by for the two disk crossings during an orbital period 
T t , and N ~ t G /T t . According to ((7j) the 2-velocity component v Z)C oc xT 1 , where 
%c — r c/L. Therefore one has to calculate the sum 9m( x c) x c with summation over all 
subsequent crossing points (odd and even) of a clump orbit with the Galactic disk. An 
important simplification in calculations follows from the fact that precession velocity 
is constant. For this reason the points of successive odd crossings are separated by 
the same angles 0. The same is also true for the successive even crossings. Using this 
simplification one can do the following transformation: 

N 



9 m {.x) 2 x 2 ~ -i- J g m (x) 2 x 2 d(f) -TT, j 9m{.x) 2 x 2< ^ dx^, (21) 



8=1 



\d>\J— I0L-/. " W dx T t 



where 



yl/2 



OlX ry.2 



x 2 \je — tp(x) — y/x 2 



(22) 



is an equation of clump trajectory in the halo. Using the above formulae one finds 



where 



AE/\E\ = (23) 
sin 7 



A(x e v) = ( 5 - jg T 9 2 m (s)ds 

1 ' 27r2(5-/5)G> o LV /2 |0|Tt2 n y/e-1>(s) -y/s 2 



3. Anisotropy of clump distribution 

A tidal heating and final destruction of clumps by the gravitational field of the Galactic 
disk depends on the inclination angle 7 of a clump orbit to the disk according to ([1]) 
and ( 1231 . This is a cause of the anisotropic clump number density decreasing during the 
lifetime of the Galaxy. A tidally induced anisotropy of clump distribution can be taken 
into account by inserting into the integral (fT5|) an additional factor e~ AE ^ E ^ = e~ ta ^ td , 
where td is an effective time of clump destruction. Due to an "anisotropy" factor 1/ sin 2 7 
in (1231) the "tidal shocking" by the Galactic disk is most effective for clumps with orbit 
planes which are near coplanar to the disk plane (that is with an inclination 7 < 1). 
Meanwhile for any intersection points there are many others orbits which pass through 
the disk plane with 7 ~ 1. For this reason the resulting anisotropy in clump distribution 
is rather small. 
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Figure 1. The normalized fractions of DM clumps in the halo P(r,a) from (|25|) 
which survive the tidal destruction by the stellar disk as a function of angle a between 
a radius- vector r and the disk polar axis. The plots are shown for radial distances from 
the Galactic center r = 3, 8.5 and 20 kpc (from the bottom to the top). The curves 
must be multiplied by factors 0.04, 0.4 and 0.9 respectively to reproduce the actual 
values. 



Let us consider a point in the halo with a radius vector r and an angle a with a 
polar axis of the disk. Only orbits with 7r/2 — a < 7 < 7r/2 go through this point. A 
survival probability for clumps can be written now in the following form 

, — 1 sin a 1 

P(x,a)= 4 2 I dp I dcos 7 / de[e-^(x)] 1/2 F(e)e- AE/lEl , (25) 
p(x) sin a J J J 



where AE/\E\ is defined by ( 123]) . This expression is derived from (Tl5l) by inserting the 
distribution over additional parameters p and 7 (by taking in mind that we use the 
isotropic halo model) and the exponential factor for clumps destruction e _A ' B// ' E '. The 
numerically calculated triple integral from (|25|) for survival probability P(r,a) is shown 
in the Fig. [TJ The annihilation anisotropy is artificially enhanced in the Fig. [T] for better 
visualization for three chosen radial distances from the Galactic center, r = 3, 8.5 and 
20 kpc respectively by using the different multiplication factors. 
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4. Annihilation anisotropy 

For the diffuse distribution of DM in the halo, the annihilation signal (e. g. gamma-ray 
or neutrino flux per unit solid angle) is proportional to 

' max (0 

'h= / pUOdx, (26) 
o 

where x = r/L and integration over r goes along the line of sight, £(C> r ) — ( r2 + r © — 
2rr cosC) 1 / 2 is the distance to the Galactic center, r max (£) = (R^ — r^sin 2 ^) 1 / 2 + 
r cosC is the distance to the external halo border, ( is an angle between the line 
of observation and the direction to the Galactic center, is a virial radius of the 
Galactic halo, r = 8.5 kpc is the distance between the Sun and Galactic center. The 
corresponding signal from annihilations in DM clumps is proportional to the quantity 

m 

''max (0 

U = »Sp J Pn (0P(Z,a)P sp (0dx, (27) 
o 

where p ~ 0.05 is a fraction of the DM mass in the form of clumps, P sp is a survival 
probability of clumps due to their tidal destructions by stars in the halo and bulge from 
[19J. The function 5* depends on the clump density profile and core radius of clump [15] 
and we use 5* ~ 14.5. Here for simplicity we do not take into account the distribution 
of DM clumps over their internal densities. As a representative example we consider 
the Earth-mass clumps M = 1O _6 M originated from 2a density peaks in the case of 
power-law index of primordial spectrum of perturbations n p — 1. The mean internal 
density of these clumps is p ~ 7 x 10~ 23 g cm -3 . The values of p and S, as well as the 
distribution of the clumps over various parameters influence the annihilation signal but 
only weakly affect the predicted anisotropy. 

In the Fig. [2] the annihilation signal calculated according to (j27|) is shown for the 
Galactic disk plane and for the orthogonal vertical plane (passing through the Galactic 
center) as function of angle £ between the observation direction and the direction to 
the Galactic center. For comparison in the Fig. [2] is also shown the signal from the 
spherically symmetric Galactic halo without the DM clumps (126]) . The later signal is 
the same in the in the Galactic disk plane and in vertical plane and therefore can be 
principally extracted from the observations. 

The difference of the signals in two orthogonal planes at the same ( can be 
considered as an anisotropy measure. Defined as 5 = (I2 — it has a maximum 

value 5 ~ 0.09 at ( ~ 39°. 

5. Discussions 

A total anisotropy of DM annihilation signal is determined in general by the Sun off- 
centering in the Galaxy and by the halo nonsphericity. An annihilation signal from 




Figure 2. The annihilation signal (f2T)) in the Galactic disk plane and in vertical 
plane as a function of the angle £ between the line of observation and the direction to 
the Galactic center. For comparison it is shown also the annihilation signal from the 
Galactic halo without DM clumps (|26|) . The values of both integrals (|27|) and (|26|) are 
multiplied by factor 10 48 . 



the Galactic center depends on the central density profile of the Galactic halo. In the 
case with a density cusp [39j, the bright source in the Galactic center is inevitable. 
Meanwhile, this cusp in the diffuse DM may be destroyed by the stellar feedback 
|40j. The small-scale DM clumps are completely destructed inside the Galactic stellar 
bulge region. The "gamma-rings" are predicted in other galaxies due to the absence of 
clumps in their centers [TT] . The unknown nonsphericity of the halo is a main source of 
anisotropy uncertainty. The value of anisotropy due to nonsphericity of the halo may 
be several times larger than one caused by the discussed in this paper effect of tidal 
destruction of DM clumps by the disk. More detailed analysis is required to separate 
these two sources of anisotropy. A nonsphericity (oblateness) of the halo due to the 
angular momentum can be easily estimated. It is natural to assume that the DM halo 
and disk have the same value of specific angular momentum (i. e. an angular momentum 
per unit mass). In this case the model of the Maclaurin spheroid for the halo gives only 
~ 0.5% difference for the halo axes. Therefore, the nonsphericity of the halo due to the 
angular momentum produces a negligible anisotropy. 

In [T7j the anisotropy with respect to disk plane in the Galaxy was pointed out. 
As it is seen from our calculations (see Fig. [2]) this anisotropy of annihilation signal 
from the DM clumps with respect to the disk is rather small, ~ 9%, but far exceeds 
the anticipated GLAST resolution, ~ 0.1%. Therefore, the discussed anisotropy may be 
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used in future detailed gamma-ray observations for discrimination of the annihilation 
signals from the DM clumps, diffuse DM in the Galactic halo and diffuse gamma-ray 
backgrounds. 

In conclusion we list some unsolved problems and unclear features of DM clump 
physics. First of all, the detailed analytical theory of small DM clump clustering is 
desirable. An effective index of the density perturbation power spectrum n — > —3 
at small-scales. This means that a gravitational clustering of small-scale structures 
proceeds very fast. As a result, the formation of small-scale DM clumps and their 
capturing by the larger ones are nearly simultaneous processes. The only approach 
developed to track this clustering is an approximate theoretical model [T5l [T9] and 
restricted numerical simulations [13]. During the past years the minimum clump mass 
is widely discussed. A value of the minimum mass has a principal significance for 
annihilation signal calculations. For this reason, the detailed calculations of minimum 
clump are requested. The central parts of DM clumps dominates in the generation of 
annihilation signal. For this reason the crucial problem is a value of the central density of 
DM clumps. To recover this value a detailed theoretical models and/or high- resolution 
numerical simulations are needed. The other problem is a survival of the central core 
of small-scale DM clump in tidal interactions. It is possible in principle that survived 
cores of DM clumps dominate in the annihilation signal. In this paper we evaluated 
only one particular part of the clump anisotropy in the halo. It would be useful also to 
clarify the halo shape and the nonsphericity of the distribution of clump orbits due to 
influence of the gravitational potential of the Galactic disk. 
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